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Abstract 

The meander problem is a combinatorial problem which provides a toy model of the 
compact folding of polymer chains. In this paper we study various questions relating to 
the enumeration of meander diagrams, using diagrammatical methods. By studying the 
problem of folding tree graphs, we derive a lower bound on the exponential behaviour of 
the number of connected meander diagrams. A different diagrammatical method, based 
on a non-commutative algebra, provides an approximate calculation of the behaviour of 
the generating functions for both meander and semi-meander diagrams. 
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1 Introduction 



The problem of enumerating meander and semi-meander diagrams has appeared in many 
diverse areas of mathematics pi and computer science [||, In the case of semi- 
meanders it can be reformulated as an enumeration of the foldings of a strip of postage 
stamps H|, |(|, which is equivalent to a model of compact foldings of a polymer. Thus the 
semi-meanders serve as a toy model of the statistical mechanics of such polymer chains. 

The generating function for meanders can be written as an Hermitian matrix model 
or as a supersymmetric matrix model R J8|. Alternatively combinatorial methods can 
be used |5|, |6|, one can represent the problem in terms of non-commuting variables || 
or using a Temper ley-Lieb algebra || |l0|, pT[| . Despite all of these different approaches 
and the initial apparent simplicity of the problem, the task of counting the meander and 
semi-meander diagrams has remained largely unsolved. 

In this paper we apply diagrammatical methods to the study of meanders. In section || 
the problem is defined and some tables of meander numbers are given. Then in section || 
we calculate bounds on the number of connected meanders, in particular the lower bound 
is improved by mapping the meanders to folded tree graphs. Section || uses a diagram- 
matic representation of a non-commutative algebra to calculate approximate generating 
functions for meanders and semi- meanders. The nature of the phase transition for semi- 
meanders is also considered in this section. Finally we give our conclusions in section [| 
of the paper. 



2 Definition of the meander numbers 
2.1 Meanders 

Suppose that we have an infinite line (or "river") on a planar surface. Then a connected 
meander of order n is defined to be a non-self-intersecting connected loop (or "road") 
which crosses the river 2n times (the crossing points are referred to as "bridges"). Two 
meanders are considered to be equivalent if it is possible to smoothly deform one meander 
into the other without changing the number of bridges during the process. The total 
number of inequivalent connected meanders of order n is denoted , and we can define 
a generating function, M(x) for connected meanders, 

oo 

M(x) = Y / M^x n . (1) 

n=l 

The first few meander numbers are MP = 1, M 2 (1) = 2 and M 3 (1) = 8. Fi gure p] shows 
some simple meander diagrams. 
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Figure 1: Some simple meander diagrams. 
The definition can be extended to include cases in which the road consists of k dis- 

(k) 

connected parts, giving a number of meanders of order n denoted by Mn . This leads us 
to define a more general generating function, 

oo n 

M(x, m) = J2J2 M n ] xnmk - ( 2 ) 

n=l k=l 

Table [l] gives M^' for small values of k and n (these numbers are taken from ||). Some 
examples of disconnected meanders are shown in figure 0. 




Figure 2: Some simple disconnected meander diagrams. 
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Table 1: Meander numbers, M« (from ref. ||). 



2.2 Semi-meanders 

Suppose that instead of having an infinite river we have only a semi-infinite river. The end 
of the river is often referred to as the "source" . In this case the road can wrap around the 
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source of the river and the corresponding diagrams are called "semi-meanders" . As before 
two semi-meanders are equivalent if one can be smoothly deformed into the other, without 
changing the number of bridges. Figure ||| shows some examples of semi-meanders. 




Figure 3: Some semi- meander diagrams. The source is marked with a dot. 

Let us denote the number of inequivalent semi-meanders with n bridges and a road 
consisting of k connected parts, as M n . Then the corresponding generating function is 

oo n 

5(c,m) = ^^MfcV. (3) 

n=l k=l 

Note that the meanders can only have an even number of bridges, whereas semi-meanders 
can have both even and odd numbers of bridges. This is the reason for the slightly different 
definitions of the generating functions and one should define x = c 2 . Table § gives the 
numbers for the simplest semi-meanders. 
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Table 2: Semi-meander numbers, M n (from ref. ||). 



3 Connected meanders 

3.1 Bounds on the number of meanders 

In this section we will derive some bounds on the number of inequivalent connected me- 
anders (from now on it should be understood that we are referring always to inequivalent 
meanders). The standard way of doing this is to consider a meander as being formed from 
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two arch diagrams (see figure || for some examples of arch diagrams). That is, if we take 
a meander diagram and cut it in two along the river, then this gives one arch diagram 
above the river and one below it (see figure ||). 




Figure 4: Two example arch diagrams. 




An arch diagram of order n consists of 2n points arranged in a line, which are connected 
up in pairs by drawing n arches above the points. The arches are not allowed to intersect 
each other and each point has exactly one arch joined to it. The number of inequivalent 
order n arch configurations is denoted A n , and the generating function for arches is defined 
as 

oo 

A(c) = J2 A nc n , A = 1. (4) 

n=0 

The empty arch diagram Aq has the value 1. A non-empty arch diagram can be decom- 
posed into two arch diagrams separated by the arch connected to the rightmost point. 
This is illustrated diagrammatically as A = 1 + A {^A \ . That is, A(c) satisfies the 
equation 

A = l + cA 2 . (5) 

Hence 

A(c) = 1 (l - VT^) (6) 

and 

An - ( 2 ") ! „ _£L. (7) 

n!(n + 1)! n 3 / 2 
In fact A n is equal to c n , the Catalan number of order n. 

Since gluing together two arch configurations of order n gives a (possibly disconnected) 
meander of order n we have 

n 

M«<^M( fe ) = V- (8) 

k=l 
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In fact, given an arch diagram on top, we can always find an arch configuration on the 
bottom that creates a connected meander (using the algorithm given in fig. 21 of ref. Q), 
which means that A n < M n l) . 

Consider the asymptotics of M n 1 ^ for a large number of bridges, 

M« ~ const. — , (9) 
where R gives the exponential behaviour and a is a configuration exponent. Then 

A n < M« < A n 2 (10) 

implies that 



2<R<4. (11) 



3.2 Improved lower bound 



The contents of the previous section are all well-known. In this section we would like to 
improve the lower bound on R. Instead of considering arch diagrams we will rewrite the 
connected meanders as tree diagrams. 

The road of a connected meander divides the plane into two parts: an inside and 
an outside. Smoothly deforming the road in such a way as to reduce the inside area 
towards zero yields, in the limit, a connected tree graph (see fig. |Bp. Given such a folded 
tree graph we can easily reconstruct the corresponding meander. Note that applying this 





Figure 6: Converting a connected meander, with the inside shaded, to a (folded) tree 
diagram. 



procedure to a disconnected meander would yield a graph that is disconnected and/or 
contains loops. The tree graphs corresponding to connected meanders can have vertices of 
arbitrary coordination number, and the mapping procedure is defined such that adjacent 
vertices on the tree occur above and below the river in an alternating fashion. There is a 
one-to-one correspondence between such tree graphs and the set of connected meanders. 
By putting bounds on the number of such trees we will improve the lower bound on R 
discussed in the previous section. 
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It is convenient to label the vertices of each tree with either an 'A' or 'B', signifying 
that they are above or below the river respectively. The left-most link in the tree will be 
marked by drawing it thicker. Each tree can then be unfolded as shown in figure |?]. 

B ABA 




Figure 7: An unfolded tree diagram corresponding to the folded tree in fig. [6|. 

In general an unfolded tree has many different ways of folding it, subject to the 
constraints that when folded the 'A' vertices are above the line and the 'B' vertices are 
below it, and that the marked link is leftmost. For example, figure can be refolded as 
in figure ||. Each unfolded tree has at least one possible folding and usually has many 
more. Thus each unfolded tree represents a whole class of connected meander diagrams. 
By counting the number of unfolded trees we will generate a lower bound on the number 



A 




B 

Figure 8: Another more regular folding of the tree diagram. 

of connected meanders. 

A rooted tree with arbitrary coordination numbers can be generated using the relation 
illustrated in fig. |9|. That is, if each link in the tree is weighted with x, then the generating 
function for trees, T(x), satisfies 

T[x) = x(l + T + T 2 + T 3 + • • •) = s/(l - T) (12) 

and hence 

T(x) = \{l~ VT^4x) • (13) 
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Figure 9: Relation which determines T{x), the generating function for rooted trees with 
arbitrary coordination numbers. 




Figure 10: Relation which determines U(x), the generating function for unfolded trees. 



Unfolded trees have a generating function U(x), which from fig. 10 is given by 



Hence 



U(x) = — = — ( 1 - VT^Ix~) - 1 = A(x) - 1, (15) 
x 2x \ > 



x 2x 

with A(x) defined through equation So we have 



U[x) = & n x n . (16) 

71=1 

However, as mentioned earlier, the number of unfolded trees (with n links) gives a lower 
bound on the number of connected meanders (with 2n bridges), that is, A n < and 
R > 2. This reproduces the result in section [O]. Note that since each link of the tree 
graph is equivalent to two bridges in the corresponding connected meander (fig. g), the 
factor for links in the tree had to be x for consistency with the definitions in section 0. 

The representation in terms of trees allows us to improve on the bound by counting 
more of the possible foldings for each unfolded tree. Suppose that the generating function 
for partially folded trees is denoted by 6(x), 

oo 

6(x) = ®nX n , (17) 
n=l 

where by partially folded trees, it is meant that only some of the possible foldings are 
included for each unfolded tree. That is, we are summing over all the trees in U(x), 
but only some of the foldings. Thus we will have A n < Q n < . Let us now define 
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0(x) = 



V) River 




Config 



Figure 11: Relation which determines G(x), the generating function for partially folded 
trees. 
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Figure 12: The relation determining T(x), which has been rewritten. 



Q(x) and hence the subset of the possible foldings that are included in the partial folding. 



First, the relation in fig. [H] is rewritten to include more foldings (see fig. 11). The rooted 
trees T(x) have been replaced by V(x), which is more folded that T and is defined below. 
The sum is over all possible ways of interleaving the arbitrary number of V^-trees. This 
relation gives 

x 



Q(x) 



1 - 2V" 



(18) 



where the x comes from the marked link. The (1 — 2V) 1 factor is due to the sum over 
all possible numbers of V-trees, with each tree being connected to one of 2 vertices ('A' 
or 'B'). 

Now consider redrawing the relation in fig. |9] as shown in fig. Then we see that 
the generating function for T only folds in one direction, but that more foldings can be 



included if we fold to both left and right. Thus for V(x) the relation in fig. 13 will be 
used instead. Figure |l4| shows a typical V-tree. 

® 





Figure 13: Relation which determines V{x) 
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Figure 14: A typical V-tree contained in V(x). 




Thus we have 

V(x) = x (l + 2V + 3V 2 + 4V 3 + • • •) = - _ X ^ )2 . (19) 
Diagrams such as that in fig. Il5| are generated when this generating function for V(x) 



is inserted in to (18). Thus Q(x) will contain a large number, but my no means all, of the 
possible folded tree diagrams. As x is increased from zero it reaches some critical value 
x c at which Q(x) is non-analytic. This is caused by the critical behaviour of V(x) due to 



4rj- diverging as x — > x c . Now, 

dV 1 



dx (1 - 4V + 3V 2 ) ' 



(20) 



so that V(x c ) = | and x c = ^. Thus 



27\ n 



4 / 

and hence R > ^ « 2.598. This improves the bound given previously. 



3.3 Further improvement to lower bound 

The lower bound can be improved further by including more foldings than exist in the 
generating function for V(x). If V(x) is replaced by W(x), whose generating function is 



10 



given by the relation in fig. 16, then an improved generating function is 

x 



9imp. (^0 



1 



where 



W(x) 



(1 



l-W> 



2W' 

(1- W) 2 
X {l-2W) 2 



(22) 



(23) 



The generating function for W{x) includes all the trees contained in V(x), but also trees 
in which there are arbitrary numbers of W-trees (emerging from the root) interleaved 
with the branches connected to the vertex. Replacing V by W in ( |l9|) and then changing 
W to W/(l — W) on the right hand side gives (|23|), This last replacement takes into 
account the fact that each branch in fig. ^ now has an arbitrary number of trees inserted 



between it and the previous branch or trunk. Figure 17 shows a typical diagram included 
in W(x). 






+ • • • 



Figure 16: Relation which determines W(x), where we are summing over arbitrary num- 
bers of W-trees inserted below each branch. 




Figure 17: Typical diagram contained in W(x). The root link of the tree is drawn thicker. 
The links which occur in a W-tree, but not a V-tree are drawn dotted. 



Equation ( ^5|) gives 



dW 

dx 



[l-Wf 



(1 - 8W + 12W 2 - 2x(W - 1)) ' 
It is the divergence of this derivative which defines the critical value of x, 

x c = ^ (71 - 17^T7) 



(24) 



(25) 
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and hence the lower bound on R is 

R > Xc -^ rs 2.970. (26) 

This is a significant improvement over the lower bound R > 2 that we had initially. It 
is worth noting that computer enumerations of the number of semi-meanders yield an 
estimated value of R m 3.50 (see reference |12|] ). A power series expansion of ®i mp .{x) 
gives 

@im P . {x)=x + 2x 2 + 8x 3 + 42x 4 + 252x 5 + 1636 x 6 + • • • , (27) 
which should be compared with the correct series for meander diagrams 

M(x) = x + 2x 2 + 8x 3 + 42x 4 + 262x 5 + 1828x 6 + • • • . (28) 

The two series agree up to and including terms of order x 4 , after which we see that ®i mp .(x) 
is, as expected, undercounting the number of meander graphs (the relevant terms in Oj mp . 
are underlined). In principle one could improve the generating functions Oi mp .(x) and 
W(x) in order to raise the lower bound further, however it appears to require increasing 
amounts of effort for diminishing returns. 

4 Diagrammatic method 

In this section, we will develop a diagrammatic method for generating meander and semi- 
meander diagrams. By truncating some of the equations, we will derive approximate 
solutions for the critical behaviour of M(x,m) and S(c,m), which can be compared with 
results from computer enumerations carried out by other authors. 

4.1 w-diagrams 

Following reference ||, let us introduce a set of non-commuting variables {u a ,ul} with 
a = 1, ■ ■ ■ ,m, which obey the relation 

u a u\ = S ab . (29) 

We will consider u a , u\ as annihilation and creation operators in some Hilbert space, with 
a vacuum such that 

uo|n) = o, (o|4 = o, <n|«) = i. (30) 

Now define G and G^ by 

G = 1 + c Gu a Gu a (31) 
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where we are summing over the repeated index a (a = 1, ■ ■ ■ ,m). Diagrammatically we 
can write the first equation as G = 1 + Gf^G^, where a line emerging from an V, that 
is written below the equation, represents u a . The fact that the two lines are joined into 
an arch indicates that we are summing over the repeated index 'a'. Since 'a' is just a 
dummy index it can be omitted. It is to be understood that each arch has a factor of c 
associated with it. Expanding this equation gives 

G = l + (^ + ^ + (^ + C\n(^ + ^(^ + (^n + £^ + ((^ + --- , (32) 

that is, G = 1 + u a u a + u a u a UbUb + u a U(,ui,u a -\ . Thus G generates arch configurations, 

however unlike A in equation (||) G keeps track of where the ends of the arches are 
positioned. The generating function G' generates arch configurations in terms of u\, 
diagrammatically one can consider these as inverted arch diagrams. So 

Gt = i + u + uu + iuj + uuu + uy + iuju + iuu; + iiyjj + • • • . (33) 

The dagger operator acting on a diagram can be thought of as reflecting the diagram, for 
example, 



it 



uy. (34) 
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Now consider the expression GG^ , this contains terms of the form u ai \ 
When this term is drawn as a diagram the inverted arches can be drawn under the upright 
arches to make the connection with meander diagrams clear. One can apply equation ( |29| ) 
repeatedly to eliminate u and v) variables. In terms of diagrams this corresponds to gluing 



together a leg on the upright arch diagram to one on the inverted diagram (see fig. 18). If 
the number of u variables is equal to the number of v) variables then we will end up with 
a (possibly disconnected) meander diagram. However, if the numbers are not equal then 
there will be spare u or v) left over. Taking the vacuum expectation value will kill these 
extra terms and hence (GG^) generates all the meander diagrams, since any meander can 
be made by gluing together two arch diagrams. Each closed loop in the meander diagram 
will give a factor of 5 aa = m. Thus the correct weightings are generated and hence 

<GG f ) = l + M{x,m). (35) 



4.2 X-diagrams 

It is convenient now to introduce a new set of variables, {X a ,X^}, defined by 

X a = ^, Xl= 2t* (36) 
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(a) 

I I ] I I I 

b, b 2 b n | b L1 b, b 2 

ii ii ii 

Figure 18: (a) Applying equation ( |29| ) to u Ql • • • u an ul ■ ■ ■ ut . (b) Applying equation 
(p9|) twice to a term in GG^. 




|X| 2 = ^A a Xt = GGt, (37) 



0=1 



So that 



|X| 2 ) = 1 + M(x,m). (38) 
The equations for X a equivalent to (|29|) and (pO|) are 

A a X fe f = -5 ab |X| 2 , (39) 
m 

x a \n) = o, (o|xj = o, (n|n> = i. (40) 

From (|l|) we have 

|X| 2 = (1 + cmX a X a ) (l + cmXlxl) (41) 
= l + cm(A a X a +AtAt) +c 2 mX a \X\ 2 Xl (42) 

where repeated indices are summed over. In a similar fashion to the u, w diagrams, this 
can be written diagrammatically, in this case as 

(tl = l + cm(n\ + \n) +c 2 m (ml, (43) 



where the ends corresponding to non-daggered operators are on the left of the vertical 
line and those for the daggered operators are on the right. Two ends are joined into an 
arch whenever the corresponding operators share a summed over index. The vertical line 
is only present in order to indicate the boundary between non-daggered and daggered 
operators. Equation (p9|) implies that 



X a X a x\x\ = I (x a X c X\xl) , (44) 

which was used above, and 

X b X a X a x\ = - (x a X a X c xf) , (45) 



m 

which can be drawn respectively as 



nin = -fftn no) 

m 
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and 

fnh = inm 



(47) 
m 



In principle one could repeatedly substitute equation ( |42| ) into itself in order to expand 
|X| 2 in terms of X a and X\. However this rapidly becomes tedious and it is necessary to 
find a simpler way of calculating |X| 2 . Let us define {A n } by 

A = l, A 1 = (\), A 2 = (ffu, ^ 3 = ^, ••• (48) 

and {Y n } by 



Yi = cm(n\ + \n), Y n = (Y n ^\. (49) 

For later use we will define the following notation, that Pi is an operator, which multi- 
plies the expression following it by X a X a (that is, Pi) on the left-hand side, whereas P 
multiplies the expressing following it by X\x\ on the right-hand side. Thus, 

Yl = cm (PI 1 + 1 Pi) = cm (p + P ) I (50) 

and 

y 2 = c (nrh + rhn) = c (n + h ) rh . (51) 

Now equation ( [43| ) can also be written as 

Ai = l + yi + aA 2 , (52) 
where a = c 2 m = xm. This equation generalizes to 

A n = A n _ 1 + Y n + aA n+1 . (53) 
Note that we wish to evaluate (Ai) since 

(Ax) = ((tl) = (|X| 2 \ = (GG^) = 1 + M(x, m). (54) 



The next step is to rewrite A\ in terms of {Y n } and a. We rewrite (p3|) as 

= a (A n+1 - fA n ) ~j( A n- fAn-i) + Y n , (55) 



af + ~ = 1, /(a) = -L f 1 - ^/T^4a") (56) 



where /(a) is given by 
and define A n by 

l n = An - (57) 
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Then from ( |55|) 

A n = fY n + afA n+1 . (58) 

This yields immediately 

A 1 = f + A 1 = f.(l + Y 1 + (af)Y 2 + (af) 2 Y 3 + •••). (59) 

The first term in (Ai) is f(a), which generates the Catalan numbers (see equation (^)). 
These correspond to the set of numbers in M(x, m), that is, the leading diagonal of 
table § 



4.3 Calculation of (Y n ) 



In order to calculate (Ai), or at least an approximation to it, we need to consider how to 
calculate (Y n ). It is worth noting that (Yi) = from (flCf). 

Now consider Y2, 



Y 2 = c (n +h)A 1 = c /(n + n ) (1 + y x + (a/)y 2 + («/) 2 y 3 + 



so that 



y 2 = (i-ca/ 2 (n +n)) [c/(n + n)(i + y 1 ) + o(y 3 ) 



(60) 



(61) 



So far all the equations have been exact, but from now on in this expression we shall 
ignore the terms Y n for which n > 3. This will make the calculation more tractable, but 
will cause us to miss out some of the meander diagrams in our generating functions. That 
is, the equation for (A\) will contain only a subset of the meanders. Now, 



Y 2 = c fY J (c«/ 2 ) n (n + n ) n+1 [1 + cm(n + n )] . I + o{y c 



n=0 

and rearranging this gives 



(62) 



Y 2 = ^- + xmf(l + xf 2 )Y / [caf 2 ) n ' (n + H I + 0(Y 



n-l 



, n+1 



(63) 



n=l 



This summation will be truncated so as to keep terms up to and including n = 2, again 
this will cause us to lose some diagrams, 



Y 2 = lh. + a;m/(l + x/ 2 ) (o + P> ) 2 + (ca/ 2 ) (n + H 



o(y 3 ) + o n +n 



(64) 



From now on the omitted terms will just be indicated by three dots, however it should be 
understood that this means terms of the form indicated in the second line of the above 
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equation. To simplify this equation, consider that, 



(n+n) 2 .l = (nn+nn).l+2(nn).l 

= (nn +nhV l+— fffo (65) 

v / m 

and similarly 

(n +n) 3 .l= (nnn + nnn).l+ — (n +n) fifrl. (66) 

m A 2 

Now from A 2 = f {A 1 +Y 2 ) + 0(Y 3 ) so that 

(n+n) 3 .l = (nnn +hhn).l+^(n +n) (A! + y 2 ) + --- 

nnn +nnnVl+-^-y 2 + — (n +n) 2 A 1 + ---, (67) 



where we can use (|59| ) to eliminate Ai. Basically one can repeat this process of substituting 
equations into (|67| ) to gradually increase the coefficient in front of Y 2 , which has the effect 
of improving the accuracy of the final result. The actual equation that is finally used is 
somewhat arbitrary, but eventually one gains an expression such as 

(n+n) 3 . 1 = (i-3x/ 2 -3x 2 / 4 ) -1 (nnn + nnn). I 



+ K Y2 + 3 J^(n + n 

mc m \ 



+ •••. (68) 



Substituting this into (|6"4|) gives, after a bit of manipulation, 



Y 2 = ^(^Y x + F.A 2 + ^(r\r\ +hn)\ + j(^c\c\ +hnn).I^J +•••, (69) 



where 



This gives us 



D = 1 - 3xf - 3x 2 / 4 (l + m(l + xf 

E = L (l _ 3 xf - 3x 2 / 4 ) , 

F = 2xf (l + xf) (l - 3xf) , 

J = x 2 m 2 / 3 c(l + x/ 2 ) . (70) 



(Y 2 ) = ^(A 2 ) + ... = a'(A 2 ) + ..., (71) 
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where we define a' = F/D. To calculate the critical behaviour it is necessary to have a 



formula for (I3). This can be derived from (^) using (49), that is, 



Y 3 = - \E.Y 2 + F.A 3 + „ ^ 



Now, 



m \ J 



+ •••• 

(72) 

(73) 



where using ( |59"|) to substitute for (Tl, and keeping only the relevant terms, we have 

if^+ilnr^ = cf((n\\n + n\(n\) + --- 

+ - = M + .... 



m 



in 



(74) 



In a similar fashion, we can extract a term proportional to I3 by the following manipula- 
tion, 



^/fSf^ + ilnr^j = 1 (fririllti + ftifnrV) = c / to I n + n I fnrY 

my y 



+ 



+ 



2^3 + 



Substituting this into equation (72) gives 



Y 3 = 1 ( £.Y 2 + F.A 3 + + ^3 J + • • • • 

D \ 2 m z / 



(75) 



(76) 



Rearranging this gives us, 



where we have defined 



H = D — -Ff 



This relation generalizes to 



l n , Jcf 3 



aA 3 + bA 2 -\ , 



F 



h- E 
b = a-. 



(77) 



(78) 



Y n = aA n + y„_i H , 



for n > 3. 



(79) 
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4.4 Meander generating function 



Now 

A n = A n _i + aA n+1 + Y n , (80) 

so that for n > 3, 

= aA n+1 + (a - + (6 + + • • • . (81) 

Let us define an operator O, which adds an outer arch to the diagram that it is operating 
on, that is, OY = X a YX\. Then since OA n = A n+ i, 

= a6 2 A 2 + (a - 1)6 A 2 + (b + l)A 2 + ■■■ (82) 

and hence multiplying by exp(yO) gives 

0= f a(4~) 2 + (o- + (6 + 1) 1 (exp(yO)A 2 \ + ---. (83) 



Keeping only the explicitly written term, the solution is 

(exp(yO)A 2 ) = Ce k y, (84) 

where 

k = ~h i 1 ~ a ~ v ' (1 _ a)2 _ 4a(6 + 1) ) • (85) 

So that (-A n ) = Ck n ~ 2 for n > 2. Using the relation (]80|) for n = 1,2 gives us the 
approximate result 

(Ai) = l + aC, C=- j± -. (86) 

1 — a — a — ak 



Expanding this using equations (|70|), ( |7q) and (pq) gives 

(Ax) = 1 + mx + (2m + 2m 2 )x 2 + (8m + 12m 2 + 5m 3 )x 3 

+(38m + 84m 2 + 56m 3 + 14m )x 4 + 0(x 5 ) + missing diagrams. (87) 

Comparing with table El we find that the approximate formula generates the correct 
coefficients up to x 4 , except for the coefficient of mx 4 , which is lower than it should be. 
Expanding the series further one sees that more and more of the coefficients are below 
the correct value, that is, as expected we are undercounting the number of diagrams as a 
result of the terms that were dropped during the calculation. 

Now the singular behaviour of (Ax) is given by the singularity in k caused by the 
vanishing of the square root. This gives a relation that determines approximately x c , the 
critical value of x, as a function of m. The quantity x c (m) has been calculated numerically 
using the equations given and the results are plotted in figure ^0. 
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4.5 Semi-meanders 



So far we have only discussed how to apply the X-diagrams to the problem of enumerating 
meanders, but in fact they can also be used to calculate semi-meanders. 

Suppose we consider a semi- meander with winding number equal to one (fig. [H]). Note 
that the winding number (denoted by w) is the minimum number of bridges that would 
need to be added if we were to extend the river to infinity on the righthand side. Me- 
anders are just semi- meanders with winding number zero and are generated by (^GG^. 
Semi-meanders with w = 1 are generated by the expression c {Gu a GG^u\G^j. That is, 
the upper part of the diagram consists of two arch configurations separated by a u a , which 
is linked to a u^ a below, that separates two inverted arch configurations. More generally a 



























Ay 





Figure 19: A typical semi-meander with winding number of one. 

semi- meander with winding number w is generated by c w (Gu ai G • • • u aw GG^u\ w 
So that the generating function for semi-meanders is 

S(c, m) = (x a Xl) + cm (x a X b xlx]) + c 2 m 2 (x a X b X c Xlx\xl) + • • • , 

which can be rewritten as 

oo 

S(c,m) = ^2 (cm) w (A w+ i) . 

Using our approximate formula for (A n ) this gives 

S(c, m) = 1 + Cm ( x H ) . 

\ 1 — mck J 

A power series expansion of this expression gives 

S(c, m) = 1 + mc + (m + m 2 )c 2 + (2m + 2m 2 + m 3 )c 3 + (4m + 6m 2 + 3m 3 + m 4 )c 4 
+(8m + 16m 2 + 11m 3 + 4m 4 + m 5 )c 5 + 0(c 6 ) + missing diagrams, (91) 

which is correct up to 0(c 5 ) except for the coefficient of mc 5 . As expected this formula 
undercounts the number of diagrams. 
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The critical behaviour for semi-meanders is caused by the vanishing of the quantity 
(1 — mck) for large m, and by the vanishing of the square root in the formula for k at 
small m (that is, m^0.6). The function x c (m) for semi-meanders has been calculated 
numerically and is displayed in figure |2(J. 



90 



80 



70 



Meander — t- 
Semi-meander — x- 




Figure 20: Plot of l/x c as a function of m for meanders and semi-meanders using the 
approximate formulae. 



4.6 Discussion 

The graph (fig. |20] ), which plots the approximate value of l/x c as a function of m for 
both meanders and semi-meanders, shows a number of interesting features. First of all it 
should be noted that the missing diagrams in the power series expansions for M(x, m) and 
S(c, m) tend to be those with lower powers of m, whereas the coefficients of the highest 
powers of m, (m n x n for meanders and m n c n for semi-meanders,) are exact. Thus one 
would expect our approximation to be exact in the m — > oo limit and least accurate in 
the m — > limit. Unfortunately, the most interesting behaviour occurs for small values 
of m. 

In reference |l2j the authors perform a computer enumeration of the semi-meander 
diagrams. They conclude that in the function S(c,m), there is a phase transition at 
m ~ 2, between a meander-like phase (with irrelevant winding number) and a phase 
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Meander — >- 
Semi-meander — x- 



4 
m 



Figure 21: Plot of R(m) = \j ^fx~ c for meanders and R{m) = \j ^fx~ c for semi-meanders, 
using the approximate formulae. 



in which winding is relevant. Figure ^ shows indications of such a phase transition 
developing at around m ~ 0.6. Below this value of m, x c for semi-meanders is equal to 
that for meanders. Above m ~ 0.6 the two curves split apart. This is caused by a change 
in the form of the critical behaviour for the semi-meanders, as indicated in the previous 
section. In the next section we will examine the behaviour of the winding number for 
semi- meanders as a function of m, and this also gives evidence of a phase transition at 
m ~ 0.6. One would expect that, as our approximation is improved by including more of 
the missing diagrams, the location of the phase transition would become clearer and be 
in closer agreement with ref. |l2|. 



Figure 21 shows our approximation to \ j -Jx~ c as a function of m. This can be directly 
compared with fig. 9 of ref. fl2"f , which calculates the same quantities, but using a different 
method. The two graphs are very similar, although it is immediately apparent that the 
undercounting of diagrams has pushed the curves on our graph lower than they should be. 
In particular, that paper derived an exact value for R(m) = 1/yfxl (for meanders) and 
R{m) similarly defined for semi- meanders. The result they gave was R(l) = -R(l) = 4, 
showing that our approximate value of R(l) = -R(l) = 3.3 is too low. Even so the general 
form of our graphs is the same as that derived in ref. [0|, lending support to the results 
in that paper. 
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In the limit m — > our approximation gives -R(O) = i?(0) = 2.6. Again this is too low 



as we have already proved in section 3.3 that R(0) > 2.970. Reference [12| on the other 
hand has i?(0) = R(0) = 3.5, which is comfortably within the known bounds. 

It should perhaps be noted that the paper [12] also showed that for m — ► oo one has 



R ~ \pm and R ~ m. This behaviour can clearly be seen in the linearity of the lines for 



meanders (in fig. |20j) and semi-meanders (in fig. |2lj) at large m. 



4.7 Winding number 

The phase transition for the semi-meander diagrams is believed to be a winding transition 
with an average winding per bridge equal to zero in one phase and non-zero in the other. 
This suggests that we should examine more closely the winding numbers of diagrams in 
the semi- meander generating function S(c,m). From (|89l) we see that 

(w) = — (cm) w w (A w+1 ) = — . (92) 

b{c,m) ^ Q ^—^^.^ 6(1 — mck) z 

Now from (||) the average number of bridges (n) is given by 

cdS 
Site 



, , c dS , , 

» = c-57- (93) 



So that 

(w) _ Cm (dsy* 

(n) (1 - mck) 2 \ dc) ' K } 

which gives a measure of the average number of windings per bridge; we expect that 
< (u>) / (n) < 1. Using a symbolic algebra program to calculate (|94|) and evaluating for 
different values of m (using our previously calculated values of x c (m)) gives us the graph 
shown in fig. This graph shows clear evidence of a phase transition occurring at about 
m ~ 0.6. As discussed in the previous section we are using an approximate generating 
function, which is known to be inaccurate for small values of m, so this result should be 
treated cautiously. However it is at least encouraging to see such clear signs of a phase 
transition, and one could hope that further effort to improve the m — > behaviour of the 
generating functions would increase the accuracy of this estimate. 



5 Conclusion 

In this paper we have shown how one can tackle the meander problem using diagrammat- 
ical techniques. In section || the problem of counting connected meanders was converted 
to one of enumerating folded trees. Using this connection it was shown that 

2.970 < R < 4, (95) 
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Figure 22: Graph of (w) j (n) plotted as a function of m, showing evidence of a phase 
transition. 



where R = -R(O) is defined by (|9|) and gives a measure of the exponential growth of 
connected meanders as the number of bridges is increased. This is an improvement over 
the usual bounds of 2 < R < 4, and is consistent with the numerical estimate 0] of 
R ~ 3.5. It may be possible to improve this lower bound, by modifying the method to 
increase the number of folded configurations that are included. Significant improvements 
seemed to be quite difficult to achieve, but the idea may be worth further study. 

In section ^| we used a non-commutative representation of the meander problem to 
calculate approximate formulae for the generating functions of both meanders and semi- 
meanders. The approximations were necessary in order to make the equations tractable 
and resulted in the loss of some diagrams, which should have been present in the generating 
functions. Nevertheless we have managed to extract graphs showing the behaviour of x c 
for the two models, which were qualitatively the same as those generated using a different 



method [12]. In particular we have shown that there is a region over which R(m), for the 
semi-meanders, is equal to R(m), for the meanders, giving an indication of the existence 
of a phase transition. Examination of the average winding number per bridge shows 
much clearer evidence of the phase transition and yields an estimated critical value of 



m c ~ 0.6. Comparison with ref. [12], which has m c ~ 2, suggests that our estimate is 
too low. It would seem to be worthwhile to attempt to improve the approximations used, 
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in order to narrow down the exact location of this phase transition. Many questions 
remain unanswered concerning the meander problem and the critical behaviour of the 
semi-meander generating function, and these must be left for future investigation. 
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